Phase-Locking of Vortex Lattices Interacting with Periodic Pinning 
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We examine Shapiro steps for vortex lattices interacting with periodic pinning arrays driven by 
AC and DC currents. The vortex flow occurs by the motion of the interstitial vortices through the 
periodic potential generated by the vortices that remain pinned at the pinning sites. Shapiro steps 
are observed for fields B^ < B < 2.25B^ with the most pronounced steps occurring for fields where 
the interstitial vortex lattice has a high degree of symmetry. The widths of the phase-locked current 
steps as a function of the magnitude of the AC driving are found to follow a Bessel function in 
agreement with theory. 

PACS numbers: 74.60.Ge, 74.60.Jg 



Vortex lattices interacting with periodic pinning ar- 
rays show a wide range of interesting commensurability 
or matching effects when the number of vortices is a mul- 
tiple or rational-multiple of the number of pinning sites. 
These pinning arrays can be created with lithographic 
techniques in which arrays of microholes or "antidots" 
and magnetic dot arrays Q can act as pinning sites. 
For small pinning sites only one vortex is trapped on a 
site as observed in transport measurements [p| , Lorentz- 
microscopy experiments || and simulations Addi- 
tional vortices sit in the areas between the pins and un- 
der the influence of an applied driving force they can flow 
between the vortices that have remained trapped at the 
pinning sites (^J^,^| The flowing interstitial vortices expe- 
rience a periodic potential caused by the repulsive inter- 
actions from the vortices at the pinning sites. The motion 
of the driven interstitial vortices is then analogous to an 
over-damped particle moving down a tilted washboard. 
With the addition of an AC driving current, interference 
effects in the form of Shapiro steps can be expected to 
occur when the frequency of the particles moving over 
the washboard matches with one of the harmonics of the 
driving frequency Qj. Recently Shapiro steps have been 
observed for driven vortices moving in samples with a pe- 
riodic array of pinning sites at twice the matching held 
B = 2B$ where B$ is the held for which there is one 
vortex per pinning site. The height of these current steps 
(range of phase-locking) strongly suggests that the vor- 
tex motion consists of the interstitial vortices moving in 
the periodic potential from the pinned vortices. Shapiro 
steps have also been observed by Martinoli et al. Q for 
vortices moving over a one dimensional periodic poten- 
tial created from a periodic thickness modulation. It has 
further been proposed that Shapiro steps can be seen for 
vortices in driven flux-transformers. 

In this work we investigate numerically and analyti- 



cally the Shapiro steps for driven vortices in thin him 
superconductors with periodic pinning arrays. The vor- 
tex lattice consists of the pinned vortices at the pinning 
sites and the sublattice of vortices that sit in the intersti- 
tial region. As a function of increasing drive we observe 
the interstitial vortices moving in one dimensional chan- 
nels between the pinning sites. With a superimposed AC 
drive we observe Shapiro steps. We find that for cer- 
tain commensurate fields, such as B — 2B<p, the system 
can be modeled as an overdamped driven pendulum with 
the associated phase locking. We find numerically that 
the widths of the steps depend on the magnitude of the 
AC driving as a Bessel function in agreement with the- 
ory. The Shapiro steps are most pronounced for highly 
symmetric interstitial vortex lattice arrangements. For 
B > 1B$ the steps vanish out due to complicated vortex 
configurations leading to nontrivial flow patterns. 

We numerically integrate the overdamped equation of 
motion for a vortex i 



f, = fT + f, 



<vp 



foe = Vj . 



(1) 



Here the total force acting on vortex i is f^. The vortex- 
vortex interaction potential is logarithmic, U v = — ln(r), 
and the force on vortex i from all the other vortices is 
fvv _ Y^j2i ViU v (rij) 0. We impose periodic bound- 
ary conditions and evaluate the periodic long-range log- 
arithmic interaction with an exact and fast converging 
sum [|12[ . The pinning is modeled as attractive parabolic 



wells with / 
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A = 1, 9 is the step function, rj^ is the location of 
pinning site k, f p is the maximum pinning force and 
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The pinning is placed 



in a rectangular array (L x ,L y ) with the ratio of the 
pinning radius r p to pinning lattice constant L y being 
r p /Ly — 0.164, close to the ratio 0.2 used in the exper- 
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iments |X(J. The pinning is placed in a 4 x 4 array and 
the initial vortex configurations are obtained by anneal- 
ing from a high temperature state where the vortices are 
liquid and cooling to T = 0. For certain parameters we 
have also considered simulations for pinning arrays up to 
10 x 10 and found only minor differences. We only con- 
sider the case for B > B^, so that the vortex motion will 
be strictly from the flow of interstitial vortices. The driv- 
ing force id represents the Lorentz force from an applied 
current. We gradually increment from zero simulating 
each DC current value for 17500 time steps (the normal- 
ized time step is dt — 0.003072) to obtain the average 
of the vortex velocities. The resulting DC force-velocity 
curve is proportional to the DC current-voltage curve. 
The AC offset is added as f a cos(u>t). We conduct a se- 
ries of simulations where the amplitude /„ is varied. In 
this work both the DC and AC driving forces will be in 
the x-direction. 

We first consider the B = 2B,p case where the intersti- 
tial vortices form a perfectly ordered square sub-lattice. 
The vortex trajectories above depinning are shown in 
the upper inset of Fig. 1 for this case. Here the inter- 
stitial vortices travel in one dimensional paths between 
the pinned vortex sub-lattice. Further, the moving inter- 
stitial vortex lattice retains the same square symmetry 
as the pinned interstitial vortex lattice. Fig. 1 shows 
typical simulation results of the voltage response V x — 
(1/N V ) Y2i=i 4 x versus an applied DC driving force at 
several different AC amplitudes for B — 2B^. The sim- 
ulation parameters are to — 1.6276 and L x = L y = 1.83. 
For zero AC driving the vortex velocities increase linearly 
with the DC driving force. With applied AC driving 
there are clear steps where the vortex velocities remain 
constant for a finite range of DC driving, indicative of 
phase-locking of the vortex motion. The widths of the 
steps depend on the magnitude of the AC drive. 

In order to demonstrate that the phase-locking of the 
interstitial vortex motion is indeed closely related to the 
well-known Shapiro steps in the AC driven pendulum 
equation j^] we first make the observation from the inset 
in Fig. 1 that the interstitial vortices are moving one- 
dimensionally along the x-direction at the symmetry line 
between the pinned vortices y = where L y is the dis- 
tance between two pinning centers along the y-direction. 
This allows us to write the equation of motion for the 
unpinned vortices as, 



± Xi -f?*( X ,y=^L) = f d + f ac 



(2) 



where we have neglected the pinning interaction and mo- 
tion in the transverse direction. We will make the addi- 
tional assumptions that unpinned vortices form a perfect 
rectangular lattice, meaning that they effectively do not 
interact due to symmetry, and that the pinned vortices 
are effectively pinned exactly to their pinning site; i.e., 
that the pinned vortices have no dynamics and form a 



perfect rectangular lattice with dimensions L x and L y . 

Under these assumptions each moving vortex obeys the 
following equation of motion [^2[: 



dt Xl L. 2^ 



sin(27rf-) 



COSh (27T^(fc + \) \ - C0S(27T-^) 



= fd + fac. (3) 



Considering only the leading term in the above sum, 
we can simplify the interaction between pinned and un- 
pinned vortices to yield the equation, 

^-x - ^sech ( TTy^ ) sin(2vr-^-) = f d + f a cos(wi) , (4) 
at L x \ L x J L x 

where we have considered only the contributions, k — 
— 1,0, and allowed for a relative error in the force of 

~ sech (^jf-^j , which is obviously small as long as L y /L x 

is not small. This equation describes the driven over- 
damped pendulum, and we can therefore apply the proce- 
dure for evaluating phase-locking ranges between a pen- 
dulum and an AC drive. 

Assuming phase-locking where the pendulum (vortex) 
moves with a frequency nu>, we insert the following ansatz 
(valid for large AC amplitudes) into the above equation, 
x(t) = XQ+nwj^-t+^j sinwi , and equate the DC compo- 
nents of the resulting expression, yielding the relationship 
between the applied AC force and the phase, 2ttxo/L x , 
for a given integer n: 



L, 



nL x ^- - —sech I tt^- ) J n ( ] sin(27r^-) = f d , (5) 
27T L x \ L X J \luL x J L x 

where J„ is the nth order Bessel function of the first kind. 
The size of the range, A/^, in fd for which the vortex 
motion may stay locked to the AC drive's nth harmonic 
is then given by the extreme values of sin(27r-p-): 



A/d = ^sech ( tt-^ 



J n 



(0) 



By conducting a series of simulations with different 
AC driving amplitudes we can compare our simulation 
results for the dependence of the step widths with those 
predicted from equation (6). In Fig. 2a we plot the widths 
of the locking ranges for the harmonics, n = 0, 1, 2, pre- 
dicted for our parameters from equation (6) (solid lines) 
and the widths of the simulated locking ranges, n = (•), 
n = 1 (o), and n = 2 (□). There is very good agreement 
between the simulation data and the predicted curves. 
We note that although equation (6) is for a single in- 
terstitial moving vortex at B — 23^ the interstitial vor- 
tex lattice is symmetric (see Fig. la) so the interstitial- 
interstitial vortex interactions cancel. We also obtain 
good agreement for the predicted widths from equation 
(6) and the simulations for the higher harmonics n > 2 
which are not shown here. We note that the agreement 
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between the simulation data and the predicted behavior 
is not expected to be exact since the force that the inter- 
stitial vortices experience from the pinned vortex lattice 
is not strictly sinusoidal. We have also tested equation 
(6) for different ratios of L x /L y by considering a rect- 
angular pinning array with L x /L y = 2. The ratio of 
the step widths for the different directions is « 52 in 
good agreement with the theoretical prediction of « 57. 
The agreement is still good when we compare the sim- 
ulated ranges of phase-locking with those predicted by 
equation (6) for the same parameters as above, but with 
L x = 2L y = 3.66. Since the vortices are forced in the 
x-direction, this is a case where the harmonic potential 
approximation made in equation (4) is not expected to 
be as good as for the square lattice case, L x = L y . Figure 
2b shows that simulations at the second matching field, 
B = 2£?0, (□ and o) show less than predicted ranges of 
locking suggesting that assumptions in the analysis are 
not well within validity. However, performing the same 
simulations, but at the matching field with one additional 
interstitial vortex (filled markers), reveals locking-ranges 
very close to what is predicted. This underlines that the 
harmonic potential assumption made in equation (4) is 
reasonable even for L x = 2L y . Closer examination of 
the dynamics at B — 2B^, (open markers) shows that 
internal modes in the moving vortex lattice are being 
excited and the assumption of cancelation of interstitial 
vortex interactions become invalid, which is responsible 
for the deviation between simulations and our prediction 
in figure 2b for B = IB^. It is, of course, important to 
emphasize that the overall features of the locking range 
is still predicted well by equation (6). 

The above analysis suggests that whenever the intersti- 
tial vortex lattice is rectangular and the interstitial vor- 
tex interactions therefore cancel, Shapiro steps should 
be observed and be well approximated by equation (6) 
when sech(nL y / L x ) is small. Square interstitial vortex 
arrangements are found at B/B^ = 2,1.5,1.25,1.0625. 
For other filling fractions the interstitial vortex lattice 
is not symmetrical and the interstitial-interstitial inter- 
actions do not cancel, leading to some deviations from 
the predicted phase-locking (the locking range is usu- 
ally smaller than predicted). This is illustrated in Fig. 
3 where we show the widths of the Shapiro steps for 
different filling fractions for a fixed AC amplitude and 
frequency. The Shapiro step widths for the different 
symmetrical vortex configurations at BjB^, = 2, 1.5 and 
1.0625 are essentially identical. For B/B^ = 1.375, and 
1.68 the steps are considerably reduced and some frac- 
tional Shapiro steps also appear. We find in general that 
for B^ < B < 2i?0, the filling fractions that produce 
square interstitial vortex lattices have the same Shapiro 
step widths as at B = 2B$. 

Interestingly for B > 2B^ we find that the step widths 
remain the same as at B = 2B,p; however, there is a 
component in V x of the steps that linearly increases with 



increasing fd- For increasing magnetic fields this linear 
increase in V x of the steps increases until B > 2.25i?0, 
when the steps can no longer be discerned. This linear 
increase suggests that only a portion of the vortices are 
phase locked. The images (not shown) from the simu- 
lations suggest that the extra vortices which have been 
added to the B = 2B^ sub-lattice cause an additional 
soliton-like motion which moves at a different speed than 
the interstitial vortices. To examine this we plot in figure 
4 the time dependent vortex velocities for two separate 
interstitial vortices along the same row at the n = 1 step 
(fd = 0.39). In Figs. 4a and 4b for B = 2B^ the sig- 
nals for the two particles are identical indicating that 
the vortices are moving in phase. In Figs. 4c and 4d we 
plot the signal from a row containing an extra vortex 
for B = 2.O625£?0. Here the same oscillation as in Figs. 
4a and 4b is seen, indicating that phase-locking is oc- 
curring; however, there is an additional lower frequency 
oscillation superimposed. The soliton like nature of this 
disturbance can be seen by noting this extra oscillation 
out of phase between the two vortices; similar to a kink 
soliton on a Frenkel-Kontorova chain |Q . 

In conclusion, we have observed Shapiro steps in the 
current-voltage characteristics of driven vortex lattices 
interacting with periodic pinning. At B — 2B^, where 
the vortex motion consists of the one dimensional flow of 
interstitial vortices between the pinned vortices, Shapiro 
steps are observed in agreement with recent experiments 
JToj l . We show that for certain filling fractions the equa- 
tion of motion for a driven interstitial vortex with a drive 
can be mapped to a driven overdamped pendulum. We 
derive the widths of the Shapiro steps as a function of rel- 
evant experimental parameters, and find excellent agree- 
ment between theory and simulations. For filling frac- 
tions where interstitial-interstitial vortex interactions be- 
come relevant the step widths are reduced. For B > 2B$ 
the steps begin to vanish due to an additional soliton like 
flow and other dynamical complexity. 
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FIG. 1. Average vortex velocity V x versus applied DC 
driving for varying applied AC driving of (from top to bot- 
tom) f a = 0.814, 0.5078, 0.22135, and for B = 2B . The 
curves have been shifted up from the curve for presentation. 
Here clear steps can be seen for finite applied AC drives. In- 
set: the typical vortex trajectories, positions (black circles) 
and pinning sites (open circles) from the simulations. Here 
the motion consists of interstitial vortices moving in one di- 
mensional channels between the vortices pinned at the pinning 
sites. Parameters are: ui = 1.6276 and L x — L v — 1.83. 
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FIG. 2. The solid lines show the predicted dependence of 
the widths of the steps in the velocity-force curves for vary- 
ing applied AC drive amplitude from equation (6) for the 
first three harmonics n = 0, n = 1, and n — 2. Markers 
represent the simulation results of the phase-locking range in 
DC current, w = 1.6276. AC amplitude, /„, is varied from 
f a = to f a « 2.93^. (a) L x = L y — 1.83 and B = 23^. 
(b) L x = 2L y — 3.66; open markers represent simulations at 
B — 23^ and filled makers represent simulations with a single 
unpinned vortex. 



4 



0.84 - 



0.42 




FIG. 3. Average vortex velocity V x versus driving force for 
constant AC driving f a = 0.5078 for varying fields (from top 
to bottom) B/B4, = 2.25, 2.125, 2.0168, 1.5, 1.375, 1.25, 1.0625 
. Each curve has been shifted for clarity. Parameters are: 
lo = 1.6276 and L x 
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FIG. 4. (a) The time dependent vortex velocity for two 
different vortices in a single row for B = 25^ in the middle 
of the n = 1 step for the same system as in Fig. 1. (b) 
The same for B = 2.06255^. Parameters are: u> = 1.6276, 
L x = L y = 1.83, f d = 0.38, and /„ = 0.442. 
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